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We present a class of spherically symmetric vacuum solutions to an asymptotically 
J> ■ safe theory of gravity containing high-derivative terms. We find quantum corrected 

Schwarzschild-(anti)-de Sitter solutions with running gravitational coupling param- 

m 

eters. The evolution of the couplings is determined by their corresponding renor- 
^^ ■ malization group flow equations. These black holes exhibit properties of a classical 

o: 

Schwarzschild solution at large length scales. At the center, the metric factor remains 
smooth but the curvature singularity, while softened by the quantum corrections, 
persists. The solutions have an outer event horizon and an inner Cauchy horizon 
which equate when the physical mass decreases to a critical value. Super-extremal 
solutions with masses below the critical value correspond to naked singularities. The 
Hawking temperature of the black hole vanishes when the physical mass reaches the 
critical value. Hence, the black holes in the asymptotically safe gravitational theory 
never completely evaporate. For appropriate values of the parameters such stable 
black hole remnants make excellent dark matter candidates. 
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I. INTRODUCTION 

One of the most challenging tasks facing theoretical physicists today is the construction 
of a consistent ultraviolet (UV) complete theory of gravity Weinberg has suggested that the 
effective quantum field description of a gravitation theory may be UV complete and non- 
perturbatively renormalizable by virtue of asymptotic safety (AS) [1]. In this scenario the 
renormalization group (RG) flows have a fixed point in the UV limit and a finite dimensional 
critical surface of trajectories approach this point at short distances. This theory has been 
extensively studied in the literature [2-9], and recent evidence suggests the UV critical 
surface is only three-dimensional even in truncations of the exact RG equations with more 
than three independent coupling parameters [10-15]. Until now, the majority of the work 
on the subject has considered significant truncations of the action, taking into account only 
the Einstein-Hilbert and occasionally cosmological constant terms. In this paper, we initiate 
the study of black hole solutions in an asymptotically safe gravity theory including higher 
derivative terms and running gravitational couplings. While we focus on the above goal, our 
methodology is easily generalized to include a broader analysis of black solutions in higher 
derivative theories. 

Our concrete starting point is a generally covariant gravitational theory with effective 
action involving a momentum cutoff p: 
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where g is the determinant of the metric tensor g^ u , R is the Ricci scalar, R^ v is the Ricci 
tensor and R^ u \ a is the Riemann tensor. The coefficients Oj (i — 0, 1, 2a, . . . ) are dimension- 
less coupling parameters and are functions of the dimension-full, UV cutoff. In particular, 
for long wavelengths we have 
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where it is useful to define the dimensionless Newton's constant Gn{p) = Gn{p)p 2 and 
dimensionless cosmological constant A(p) = A(p)p~ 2 . The couplings satisfy the following 



RG equations, 
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The conditions for asymptotic safety require that all the beta functions vanish /3i = 
when the coupling parameters g{ approach a fixed point g*. If g* = 0, the fixed point is 
Gaussian; if g* ^ 0, the fixed point is Non-Gaussian (NG). For the NG fixed point, all the 
coupling parameters are fixed, the cutoff p becomes irrelevant as p — > oo, and the theory 
is adequately described by a finite number of higher order counter-terms included in the 
effective action. Near the fixed point we may Taylor expand the beta functions in a matrix 
form 

P l [g] = Y J B l3 {g 3 -g*) , (4) 
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where the elements of the matrix are defined by B^- = vj' — at the fixed point. Solving the 
RG equations (3) in the neighborhood of the fixed point we find 

w.) • (5) 

where e n and v n are the suitably normalized eigenvectors and corresponding eigenvalues of 
the matrix Bij. Since B is a general real matrix with symmetry determined by a particular 
gravity model, its eigenvalues can be either real or in pairs of complex conjugates. As a 
consequence, the dimensionality of the ultraviolet critical surface is equal to the number 
of eigenvalues of the matrix B, of which the real parts take negative values. The above 
solution involves an arbitrary mass scale M*. By requiring the largest eigenvector of order 
unity, M* is typically identified with the energy scale at which the coupling parameters are 
just beginning to approach the fixed point. 

It is well known that the quantization procedure for General Relativity leads to a non- 
renormalizable quantum field theory, where an infinite number of terms have to be fixed 
to renormalize the standard perturbation theory. Including higher derivative terms in the 
action (such as we have done in (1)) introduces higher derivative propagators that soften 
the divergences of the perturbative quantization and can result in a perturbatively renor- 
malizable theory [16]. Unfortunately, the higher derivative terms required, introduce new 
massive spin-two degrees of freedom with wrong-sign kinetic term (i.e. negative norm states, 



or ghosts) [17]. For example, if we quantize this system in a canonical form the problem 
with ghosts arises since the metric field and its high- derivative term are regarded as two 
independent variables. These instabilities generically render the theory non-unitary and un- 
stable. However, as shown in [18], if the quantum theory is based on a path integral which is 
evaluated in Euclidean space and then Wick rotated to Lorentzian space, the path integral 
can yield a theory of quantum gravity without a negative norm state. A specific application 
of this approach in inflationary cosmology provides an interesting interpretation for metric 
perturbations [19]. Surprizingly, evidence is mounting that these catastrophic ghosts may 
be eliminated in higher derivative theories when the UV limit of gravity is controlled by a 
NG fixed point as in the AS gravity scenario [12, 20-23]. 

Because the asymptotically safe gravity theory is relevant for the physics of high en- 
ergy scales and short distances, one is naturally lead to consider its application to early 
universe cosmology. For example, Weinberg has recently argued for the existence of infla- 
tionary solutions complete with graceful exit in the context of the theory [15]. Another 
clear application of AS gravity is to black hole physics; in particular, to understand how the 
theory modifies the conventional Schwarzschild black hole solution by taking into account 
the quantum corrections naturally incorporated into the model. Some early attempts are 
presented in [24-27] . Their key assumption was that the leading order quantum corrections 
to the black hole spacetime are captured via a simple running of the Newtonian constant Gn 
determined by the renormalisation group equation for gravity; however, the consistency of 
the resultant modified metric and modified effective action was not thoroughly investigated. 
For simplicity, the analysis was limited to a severe truncation of the action (1), including 
only the gi term. By truncating to the Einstein-Hilbert term (as in any f(R) truncation) 
the authors eliminated the four-derivative propagator for the helicity two states which may 
drastically alter the AS theory and its solutions. In this paper, we develope an effective 
method of finding vacuum solutions to Einstein's equation derived from the AS gravity with 
higher derivative terms. We focus on spherically symmetric spacetimes and present an exact 
form of a Schwarzschild- (anti)-de-Sitter (SAdS) solution. The thermodynamical properties 
of these black hole solutions are briefly discussed. 

The paper is organized as follows: in §11, we derive the equation of motion for the AS 
high- derivative gravity and study a generic vacuum solution to this theory which preserves 
static spherical symmetry. We find that, to leading order, this solution is exactly SAdS. 



Due to the quantum nature of renormalization, the coupling coefficients vary with respect 
to the running of energy scale and it is necessary to investigate the RG flow of these coupling 
coefficients in accordance with the asymptotically safe scenario (see, §111). In §IV, we analyze 
the black hole solutions to the theory. We determine the scale identification between the 
momentum cutoff and the radial coordinate and subsequently, we study the behavior of the 
background geometry in the UV and infrared (IR) limits. We find that the metric factor is 
nonsingular for all values of the radial coordinate but the curvature diverges at the origin. 
Moreover, there exists a critical mass for the background geometry. When the black hole 
mass is heavier than the critical value the black hole has two horizons, but this solution 
vanishes when the mass is subcritical, revealing a naked singularity. Numerical calculations 
are performed to confirm our analytic results. The black hole thermodynamics are briefly 
discussed and we present our concluding remarks in §V. 



II. STATIC SPHERICALLY SYMMETRIC VACUUM SOLUTIONS 

Black hole physics provides a window into the quantum nature of gravity. Black holes 
possess many remarkable properties, for example, the associated thermodynamics [28] and 
holographic properties. The simplest black hole solution in four dimensions corresponds to a 
Schwarzschild spacetime [29, 30]. As is well known, the geometry is divided into two causally 
independent regions by an event horizon located at the Schwarzschild radius r s = 2G^M l . 
Given the possibility of an asymptotically safe gravitational theory, it is logical to re-examine 
these intriguing properties of black hole physics when quantum corrections are incorporated 
into the gravitational action. We now obtain the field equations from the action (1) and 
search for static, spherically symmetric solutions to the theory. 

The generalized Einstein field equations are obtained by varying the action (1) with 
respect to the metric tensor g^: 

where G^ is the generalized Einstein tensor and T^ v is the energy-momentum tensor of the 
background matter and is vanishing due to our vacuum ansatz. The generalized Bianchi 



1 Here and throughout we use Natural units with c = H = ks = 1. 
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identity is preserved automatically due to the general covariance of the scalar-type action. 
The form of G pu appearing in the above equation can be expressed as, 

&v = 1 ^ Li gQ + p 2 giR + g2aR 2 + g2b R ffp R°P + g 2c R apXK R° pX A - p 2 9l R pv 
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In the above, we have defined the d'Alembertian operator □ = V^V^, where V M is the 
covariant derivative. The trace of the Einstein tensor gives the useful quantity, 

G = g^GT 

= 2p 4 g + p 2 9l R - 2(3g 2a + g 2b + g 2c )DR , (8) 

which vanishes in our vacuum investigation, T pv = 0. In the above derivation we have 
applied the Bianchi identity: V a {R ap - \g ap R) = 0. 

We assume a static, spherically symmetric metric ansatz, 

ds i = -A(r)dt 2 + ^- + r 2 dtf 2 . (9) 

Substituting the metric (9) into the action (1) and neglecting the time integral, up to order 
O(p ) gives, 

V P [A,B] = JdrA-K^^^ + g^R + g^R 2 

+g2 b R,uR^ + g2 C R,u Pa R^ pa ) , (10) 



where the evaluated curvature invariants in terms of the metric (9) are included in the 
Appendix. To obtain the vacuum solutions we must solve: 

^ = ^ = 0. (11) 

SA SB y ' 

The above calculation is greatly simplified by choosing the Schwarzschild gauge B(r) = 
N(r)A(r) and setting N — 1. This relation is ensured by the Cauchy theorem for a Rie- 
mannian geometry with unique boundary. 
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Varying the action with respect to the function A(r) gives: 

4(3(? 2a + g 2b + g 2c ){-2 + 2A- r 2 A") 

+ (2fe + g 2b + 2g 2c )r 3 (4A^ + rA&) 
= . (12) 

Solving the above equations of motion yields, 

A(r) = 1 + Cl r ni + c 2 r n2 + — + c 4 r 2 , (13) 

where the exponents rij (j = 1, 2) are given in terms of the couplings g 2 k (k = a, b, c), 
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This solution involves four coefficients c\, c 2 , C3 and C4 which must be determined by other 
constraint equations, boundary conditions and the consistency relation with infrared limit 
solution. 

To determine the coefficients appearing in the solution (13), we impose the constraint 
equation G — 0, consistent with our vacuum assumption. A natural choice is to set c\ = 
c 2 = 0, giving the SAdS solution to the high- derivative gravity theory: 

4 . . 1 2G P M r 2 
A(r) = 1 - -^- - - , (14) 

where G p is a gravitational coupling which is varying along with the running of the cutoff 
scale. The parameter M is an integration constant we identify with the physical mass of 
the black hole, and l p is the radius of the asymptotic (A)dS space which is also a cutoff 
dependent function, expressed as 

,2 %i 
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The solution is consistent with the usual vacuum solution from GR when the energy scale 
flows to the IR limit, but differs significantly at high energy scales, as we shall see. We 
require a scale identification between the coefficients of the solution and the cutoff scale p 
which is dependent on the details of the AS gravity theory. Therefore, in the subsequent 
section we apply the AS scenario of the high-derivative gravity action. 



III. ASYMPTOTICALLY SAFE HIGH-DERIVATIVE GRAVITY 

The functional gravitational RG equation is based on a momentum cutoff for the propa- 
gating degrees of freedom and captures the nonperturbative information about the gravita- 
tional theory. The RG equation is of the form [4]: 

i£r P r ' = H smr > +R{p) )^r P n{ph (16) 

where R is an appropriately defined momentum cutoff at the scale p and is usually deter- 
mined by the so-called optimized cutoff process [31, 32]. In the above formula, we suppose 
the gauge fixing terms have already been included, although they are irrelevant for our 
present consideration. Additionally, the trace in the RG equation stands for a sum over 
spacetime indices and a loop integration. Our philosophy then is to effectively integrate out 
the high momentum fluctuations with momentum larger than the cutoff p, and incorporate 
them via the modified dynamics for the fluctuations having momentum less than p. 
The high-derivative terms in the effective action (1) can be regrouped as follows, 

(17) 
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E = R^RTi" - AR^R^ + R 2 , (19) 

is the integrand of the Gauss-Bonnet term which is topological in 4-dimensional spacetime. 
The re-expressed high- derivative terms in Eq. (17) are equivalent to the terms appearing in 
the original action (1), under the following identifications, 



92a = 


1 9 oo 
~6A + A + 3A 


92b = 


1 49 
A _ T' 


92c = 


1 9 

~2A + A ' 



(20) 

Note that a standard derivation of dimensional regularization with d = 4 — e [33, 34] 
indicates that the beta functions of the dimensionless coefficients introduced in (1) are 
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sensitive to the dimensional corrections for the case, e < 0. Therefore, we restrict our 
attention to the case of e > 0. Explicitly in the nontrivial limit of e — > 0, the beta functions 
for the dimensionless coefficients of the high-derivative gravitational terms are given by [33- 
35], 

* = -(k-^* • (21) 

To perform the stability analysis, we note that the coefficient A has the familiar logarithmic 
form which approaches asymptotic freedom, 

KP) = l + ^Aolnp/M p ' (24) 

where Ao is a fixed value of the coefficient A at the Planck scale. Since the above form is 
logarithmic, we find A ~ Ao in a wide range around the Planck scale. We will apply this 
approximation in deriving the coordinate dependent cutoff scale later. Moreover, the other 
two parameters 9 and oj, approach a group of fixed points in the UV limit, among which the 
stable ones take the values 

9* ~ 0.327 , tu* ~ -0.0228 , (25) 

where the superscript "*" denotes the parameter value at the NG fixed point as introduced 
in the beginning of this paper. 

By solving the beta functions for the gravitational coupling and cosmo logical constant, 
we observe a Gaussian fixed point in the IR limit and a NG fixed point in the UV limit. 
The central result is 

(Am + VP 2 G N )(1 + & 2 G N ) 
8np 4 G N 
± + & 2 G N 
167rp 2 Gjv 

where G^ and Km are the values of the gravitational coupling and the cosmological constant 
in the IR limit which are determined by astronomical observations. In the UV limit, the 
coefficients of the Einstein-Hilbert part of the gravity action will flow to a NG fixed point 
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with g —7- Pq and g\ — > g\. To obtain the NG fixed point values of the remaining variables 
in the high-derivative terms we solve the flow equation explicitly 

9*o * -^ * -6.331 x 10- 3 , (28) 

g* ~ — ~ 1.432 x 10" 2 , (29) 

yi 16?r V ; 

by a numerical computation. Notice that the running gravitational coupling G p is related 

to the coefficient g\ by, 

Gp = 16^ = 1 + ^G N ■ (30) 

From the analysis of g±, we see that G p coincides with the Newtonian constant Gn at low 
energy scales but decreases rapidly as the momentum cutoff p, which implies a weakening 
of gravity at high energy scales. We will see that this weakening of the gravitational force 
in the UV leads to a softening of the singular behavior of the black hole solution near the 
origin. 

IV. BLACK HOLES & ASYMPTOTICALLY SAFE GRAVITY 

From the above analysis, we find that the spherically symmetric vacuum solutions to 
the high-derivative gravity theory flow to the classical SAdS geometry at low energy scales. 
However, in order to implement quantum corrections to the running coefficients appearing 
in the classical geometry, we must determine the relationship between the momentum cutoff 
p and the radial coordinate r. 

A. Relevant scale 

Recall the quantity G, introduced in Eq. (7), vanishes for the vacuum state, T^ u = 0. 
Inserting the solution (14) into Eq. (7) and making use of the redefinition (20), we could, 
in principle, solve for the momentum cutoff pas a function of the radial coordinate r by 
requiring G — 0, which gives 

G = 2Ao + P 2 9iR ~ ^-OR = . (31) 

Note that the above equation is a fourth-order differential equation and will yield a class of 
solutions for the momentum cutoff p(r). Consequently, it is necessary to impose the relevant 
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constraints in order to determine the physical solution. We now focus our attention to the 
black hole solutions at high energy scales and at small radial distances compared with the 
Planck scale, l p . 

As a consequence of requiring G = 0, we derive a scale identification between the mo- 
mentum scale p and the coordinate r of the form, 

p(r) ~ 2.663 [ TV ^] r _ 4 , (32) 

when the energy scale is as high as the AS scale. We note that this relation is different 
from the result p ~ r ~ 3 / 2 appearing in [27] through a UV matching. The discrepancy 
is due to the action truncation considered in [27], which consisted of only the Einstein- 
Hilbert term and a vanishing cosmological constant. In the present analysis, we include the 
higher derivative terms and nonzero cosmological constant in addition to the EH term and 
running G^. Consequently, we consider the effective action in the low energy limit. When 
the momentum cutoff flows to the IR regime, the high-derivative terms are suppressed 
automatically and thus becomes negligible. In this case, p ~ 1/r which is consistent with 
the IR matching of [27]. The asymptotically safe SAdS geometry is obtained by replacing 
the classical gravitational constant Gn and the (A)dS radius with a RG improved G p (r) and 
l p (r), respectively. 

Note that the dS radius l p has significant cosmological implications and may be related to 
the nearly exponential acceleration of our universe at early times, as well as, the cosmological 
acceleration occurring today. The running of l p can drive an early inflationary period, 
complete with a successful graceful exit [15]. For late time cosmological acceleration, we 
may consider the possibility that the value of l p in the IR limit is of order of the size of 
our universe in accordance with the latest cosmological observations. Such cosmological 
applications are beyond the scope of the present analysis. 

B. Black hole solutions 

Let us focus our interests on the local quantities of the quantum corrected black hole 
solutions (especially near the horizon and the origin) and leave the investigation of dS 
radius to future study. One can finely tune a very small value of Ao in order to let the dS 
radius be much larger than the black hole horizon. In the following we study this geometry 
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along with the momentum cutoff p flowing from the UV regime to the IR regime. 



1. UV limit 

As the radial coordinate approaches the origin, the energy scale of the cutoff could be 
very high, and consequently the coupling parameters may already have arrived at their NG 
fixed point. Correspondingly, the gravitational coupling is expressed as 

Gp(r) = 16np2(r) gi ~ £pHr) ' (33) 

when p ^> M p . In this case, the coefficient A approaches zero, leading to infinitely large #2- 
According to Eq. (15), the radius of the asymptotical dS space diverges. Therefore, we can 
neglect the last term of the solution (14) near the origin. Inserting Eqs. (32) and (33) into 
(14), we derive the approximate form of the metric factor in the UV limit: 

A uv (r)^l-^-\X \l(Mr)^ . (34) 

The above result indicates that the metric factor A(r), unlike its GR counterpart, is no 
longer singular inside the horizon of the black hole in AS gravity with leading order high- 
derivative terms. Naively, it may appear that the geometry is glued to an asymptotical 
Minkowski spacetime near the origin. However, substituting the solution into the expression 
for the Ricci scalar, we obtain 

which is singular at the origin of the spacetime. By applying the horizon condition Ajjy(r) = 
0, we obtain an approximate solution 

0-6717T 2 

ruv ~ 1T7^7 ' (36) 

at high energy scales, indicating that the curvature singularity is hidden behind the UV 

horizon. The other familiar invariants are altered as well and reflect the softening of the 

singular behavior in the asymptotically safe gravity theory: 

n wx, „ 31640625 VApM 

" vXa ~ 4194304 vr 2 r3 ' [6<) 

r r ^ ~ H71875 7a m 

G ^ AffG " 4194304 vr 2 r* ' [ ^ } 

„ 59765625 V\ M 
^ ~ 8388608 vr2 r* ' l } 
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These should be compared with their hard singularity partners in General Relativity, for 
example, the Kretschmann scalar (37) in GR diverges at the origin as ~ 1/r 6 . 

2. IR limit 

Along with the increasing of the radial coordinate the momentum cutoff could drop to 
the regime below the Planck scale as p ~ 1/r. In this case, 

G p (r) ~ G N ( 1 - &L j , (40) 

where £ is a constant which deviates from the coefficient £ by a factor of 0(1). Substituting 
this gravitational coupling parameter into the black hole solution yields 

, , , , 2G N M / £G N \ 
A IR (r) ~ 1 - — — ^1 - 5_ 1 . (41) 

By solving the equation of horizon condition, we obtain another real solution which 
corresponds to the event horizon at low energy scales. Expanding to the leading order in 
the coefficient £, we obtain the familiar Schwarzschild horizon value, 

r IR ~ 2G N M - JL + (f ) . (42) 

Hence, the horizon location of the AS black hole is in precise agreement with the usual form 
for the Schwarzschild radius r s = 2G^M in the limit £ — > 0. 

By comparing the solution for the event horizon in the UV limit (36) and the horizon in 
the IR limit (42), we see there is a critical value for the mass parameter of the black hole: 



i 



M c ~ |Ao|"GV • (43) 

When M > M c , the black hole in AS high-derivative gravity has two horizons. There is 
an inner Cauchy horizon as well as an outer event horizon corresponding to that of the 
ordinary black hole solution in GR. At the extremal value M = M c , these two horizons 
coincide; beyond the extremal value (M < M c ) the black hole vanishes, corresponding to 
a naked singularity. The double horizon feature of our solution is in agreement with the 
results obtained in the AS gravity with Einstein truncation analyzed in [24]. 

Here we make further comparison between our results and those obtained in Refs. [24-27]. 
First, due to the presence of higher derivative terms, the form of the scale identification (32) 
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is altered compared with the result obtained in the case of the pure Einstein truncation [27]. 
Furthermore, the black hole solutions obtained with Einstein truncation are smooth and 
nonsingular, but become singular if higher derivative terms are taken into account. The 
conditions identifying the critical mass are the same in the two scenarios, i.e., both the 
metric factor and its derivative with respect to the mass are vanishing. 

3. Numerical analysis 

In the above analytic study we made a number of approximations in order to obtain 
the black hole solutions in the AS high-derivative gravity theory. To better understand 
the analytic results, we now turn to a quantitative numerical analysis. Our results for 
the coordinate-dependent momentum cutoff p(r) and the quantum corrected Schwarzschild 
metric factor A(r) are presented in Figs. 1 and 2, respectively. 
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FIG. 1: Running of the momentum cutoff p as a function of the radial coordinate r. Here we take 
Ao = Gn = 1- 



Fig. 1, shows the slope of the running momentum cutoff curve in the IR and UV regimes. 
In the IR regime, which corresponds to a large length scale, p scales approximately as 1/r. 
However, in the high energy limit, the relation between the momentum cutoff and the radial 
coordinate becomes p ex r -3 / 4 which is in agreement with the analytic result obtained in Eq. 
(32). 

In Fig. 2, we plot the goo metric component of the spherically symmetric vacuum solutions 
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FIG. 2: The metric factors A of the spherically symmetric vacuum solutions as functions of the 
radial coordinate r in both Einstein gravity and the AS high-derivative gravity. In the numerical 
computation, we take Ao = G/v = 1- We consider M = 1, M = 0.403 and M = 0.1, which 
correspond to a RG running black hole, a critical black hole, and a naked singularity, respectively. 

of AS high-derivative gravity for the values of the mass parameter: M — 1, M — 0.403 and 
M = 0.1. In order to illustrate the significance of the quantum corrections to the vacuum 
solutions, we compare them with the metric factor of a black hole in classical Einstein gravity 
as denoted by the black solid curve. The red dash curve corresponds to the metric factor of 
a double-horizon black hole with M — 1 in AS high-derivative gravity. When the physical 
mass decreases to the critical value (M c = 0.403 in this example), we obtain an extremal 
black hole and the outer horizon coincides with the inner horizon. This case is plotted by 
the green dash-dot curve. Moreover, the blue short dash curve shows that there are no black 
hole solutions (i.e. there is a naked singularity) for subcritical values of the mass M < M c . 
In all of the above solutions the metric factors of the background geometries are smooth and 
nonsingular, but the curvature is divergent at the origin. 

C. Thermodynamics of the IR horizon 



From the above analysis, we have shown how the normal Schwarzschild geometry is 
corrected by nonperturbative effects in quantum gravity. To further probe the quantum 
gravity effects we now study the process of black hole evaporation via Hawking radiation [28]. 
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We apply the Euclidean path integral approach to determine the temperature and the specific 
heat [36], and employ the method of complex path analysis to study the emission rate of 
the Hawking radiation [37, 38]. We now investigate the thermodynamical features of the 
quantum corrected black hole around the IR horizon in the super-critical solution with 
M> M c . 

There has been considerable interest in the study of thermodynamics of quantum cor- 
rected black holes. For example, the Hawking radiation of a two dimensional nonsingular 
black hole [39] in the context of dilaton gravity [40, 41]. Similar to the derivation performed 
in [39], we compute the Bekenstein-Hawking temperature 

1 



T = T ^ir IR ) 
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8nG N M 



AG N M 2 



(44) 



around the IR horizon. In comparison with the standard result, the temperature of the 
quantum corrected black hole is smaller than that of the classical analog (see Fig. 3). Stefan's 
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FIG. 3: Approximate Hawking temperature as a function of black hole mass (see Eq. (44)). The 
AS high-derivative solution is the solid (red) curve and is compared with the standard Hawking 
temperature of General Relativity T _1 = 87rM, dashed (black) curve, which diverges for small M. 
In the numerical computation, we take £ = .65 and Gm = 1- The temperature reaches a maximum 
at M = \/3£/4Gn corresponding to T~^ x = 6tt\/3Gn^, and reaches zero at the critical value 
M c = \IZ/AG N ~ .403. 
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law gives the black hole radiation power law derivation: 

_ dM 7 . 



" L 3i (45) 



256tt 3 G 2 v M 2 \ 2G N M 

where a = 7r 2 /60 is a Stefan-Boltzmann constant. Hence, the evaporation of the quantum 
corrected black hole can only take place when its physical mass is heavier than the critical 
value. Consider Ao ~ 0{1), this critical value coincides with the critical mass M c as intro- 
duced in the previous subsection. Once the physical mass is smaller than M c , the absorption 
process will dominate over the emission process of the black hole radiation. Therefore, a crit- 
ical mass black hole with M ~ M c is expected to be the most stable. Indeed, from Fig. (3) 
we see that the extremal black hole ceases to decay because the temperature vanishes for 
the critical mass. 

From the radiation power law obtained above, it is possible to estimate the evaporation 
time of the black hole taking into account the quantum gravity corrections 

M T S 



\V\ I 3|_ 

- 1 - 2G N M 



(46) 



where Ts is the evaporation time of a classical black hole. The quantum corrected evaporation 
time is infinite when the physical mass of the black hole decreases to the scale around the 
critical mass M c , in accordance with our expectations. 

In the above analysis we studied the black hole thermodynamics in the framework of the 
AS high-derivative gravity theory using perturbative methods with £ < 1. However, we do 
not expect these qualitative results to be drastically altered in the nonperturbative regime. 
This is due to the fact that there is always a critical mass for the black hole solution where its 
IR horizon coincides with the UV horizon. This corresponds to an extremal black hole with 
a vanishing temperature (since T oc A' = 0). As a consequence, the black hole evaporation 
time will be infinite as discussed above. 

V. CONCLUSIONS 

In this paper we have initiated the study of static spherically symmetric vacuum solutions 
to the theory of asymptotically safe gravity with high-derivative corrections. We find that 
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a generic solution corresponds to a quantum corrected SAdS spacetime due to the RG flows 
of the gravitational coupling parameters. Under certain simplifying assumptions, we obtain 
a black hole solution with smooth metric factor and a curvature singularity at the center 
of the geometry. The singular behavior is mildly softened in the AS high derivative theory 
compared to the corresponding GR black hole. The quantum gravity corrected black hole 
solution generically possesses two horizons, which correspond when the physical mass of 
the black hole decreases to a critical value. The temperature of the black hole is exactly 
zero when its physical mass reaches the critical value. We find that the temperature of the 
quantum corrected black hole is, in general, lower than that of a classical black hole and 
stable critical mass remnants are natural final states of the Hawking evaporation process. 

These results are related to many other interesting issues which deserve future study. 
For example, a collection of critical black holes may serve as an alternative candidate for 
dark matter if their masses are as low as the TeV scale. From Eq. (43), we find that this 
is a possibility if Ao ~ 10 64 which is also compatible with the theoretical requirement of 
suppressing gravitational high-derivative terms. The passive thermodynamic properties of 
the AS black hole with a critical mass may leave significant signals on the matter power 
spectrum at small length scales. In addition, if sufficient numbers of critical mass holes 
were produced at the end of inflation, they may leave potentially measurable signals in 
the polarization of the cosmic microwave background radiation on small angular scales. A 
potential concern is that if the critical mass of a black hole is much lower than the Planck 
scale Minkowski spacetime could become unstable due to the perturbation modes which 
increase the mass of the black hole. It is then possible that there would be too many 
primordial black hole remnants produced in a static background. We note that this problem 
may be alleviated in realistic cosmological models due to the expansion of the universe. As 
the expansion of the background spacetime dominates over the production rate of black 
holes, the density of primordial black holes may be significantly diluted after a sufficiently 
long period. We leave the study of this intriguing possibility for future work. 
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Appendix: Curvature invariants 

The curvature invariants appearing in the action (1) evaluated in terms of the line element 

(9): 

AA 2 - AA 2 B - ArABA' + r 2 BA' 2 - ArA 2 B' - r 2 AA'B' - 2r 2 ABA" , A . 

R = 2^A 2 (A1) 



R 2 = — j_ ( - (r 2 BA 12 ) + rA (2rBA" + A' (rB' + AB)) 

+ AA 2 (rB' + B-l)) 2 (A.2) 



R ^ RIW = —- l {r 4 B 2 A' 4 + AA 4 {A{-l + B) 2 + Ar{-l + B)B' + Zr 2 B' 2 ) 

+ ArA\A'(A(-l + B)B + 2rBB' + r 2 B' 2 ) + 2r 2 J B J B'A") 

+ r 2 A 2 (A' 2 (12B 2 + r 2 5' 2 ) + 4r£A'(25 + rB')A" + 4r 2 B 2 A" 2 ) 

- 2r 3 ABA' 2 (rA'B' + 2B(A' + rA"))) (A.3) 



R^R" uXa = ^-^(r 4 B 2 A> 4 + 8A 4 (2(-l + B) 2 + r 2 B' 2 ) 

+ -2r 4 ABA' 2 (A'B f + 25A") + A 2 (r 2 A /2 (8B 2 + r 2 5' 2 ) 

+ 4r 4 5A / 5'A ,/ + 4r 4 5 2 A" 2 )) (A.4) 



